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We briefly review recent developments of hydrodynamics, its gravitational description and 
relevance to relativistic heavy ion collisions. We discuss the basics of hydrodynamics, the 
fluid/gravity correspondence, triangle anomalies and chiral effects, turbulence and its universal 
structure. 

1 Introduction 

According to the Holographic Principle^, the degrees of freedom of a quantum theory of gravity 
in a volume of space V are encoded on its boundary A. Thus, the Holographic Principle relates 
gravity in d + 1 space-time dimensions to field theorv without gravity in one lower dimension. 
In the framework of the AdS/CFT correspondence'^lSI^ the quantum theory of gravity is string 
theory on asymptotically Anti-de-Sitter (AdS) space, and the theory on the boundary is a 
conformal field theory (CFT). The AdS/CFT correspondence has been generalized, and one has 
by now a large number of examples of non-AdS/non-CFT relations. The non-CFT theories are 
QCD-like gauge field theories. In the regime of a large field theory coupling A ^ 1 the gravity 
description is weakly coupled, that is the curved geometry has a small curvature. On the other 
hand, when A <C 1 the field theory is weakly coupled, while the gravity description is strongly 
coupled. 

A typical curved metric that arises in this setup takes the warped form 

ds^ = dr^ + a{rf{-df + dx'dx') . (1) 

Here {t,x^) are the field theory space-time coordinates, and the radial coordinate r is interpreted 
as the field theory energy (RG) scale. The warp factor o(r) encodes information about the nature 
of the field theory. When a(r) = exp(— r), the metric ([T]) is that of AdS space. A particulary 
interesting regime of the field theory is the hydrodynamic approximation, where the AdS/CFT 
correspondence relates the field theory hydrodynamics to perturbations of black hole (brane) 
gravitational backgrounds. 

An important experimental framework to which the correspondence has been applied is the 
description of the QCD plasma produced at RHIC and LHC. This plasma seems to exhibit a 
strong coupling dynamics, as{TpiHic) ~ 0{1). Nonperturbative methods that can be used to 
study real time dynamics are largely unavailable, while lattice QCD methods are inherently 
Euclidean. The AdS/CFT correspondence provides a real-time nonperturbative framework, 
and one typically uses the strong coupling properties of a CFT plasma as a reference point for 
describing the strongly coupled QCD plasma. 



A particulary studied quantity is the ratio of the shear viscosity rj to the entropy density 
s, which for a generic strongly coupled gauge field theory is low. Indeed, such a low ratio is a 
generic property of the gravitational description^. 

The hydrodynamic simulations at low shear viscosity to entropy ratio are consistent with 
RHIC dataE^ see figure [1]) , where the elliptic flow parameter is the second Fourier coefficient 
V2 = (Cos(2</>)) of the azimuthal momentum distribution dN/d(j) 

dN 

— ~ 1 + 2v2Cos{24)) . (2) 
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Figure 1: Comparison between hydrodynamic simulations and experimental values of the elliptic flow parameter 

(Luzum,Romatschke:2008). 



2 Basics of Hydrodynamics 

In the hydrodynamic regime the system is in local thermal equilibrium enforced by frequent 
collisions between the particles. Thus, it is characterized by a short mean free path (correlation 
length), much smaller than the characteristic scale of variations of the macroscopic fields. Since 
it is dominated by collisions, an appropriate description is that of a collective fluid-type flow 
rather than a particle (kinetic theory) one. 

The Hydrodynamic Modes: The effective degrees of freedom in the hydrodynamic regime 
are charge densities p{x,t), which are non-equilibrium thermal averages {p{x,t))non—eq- 

The 

hydrodynamics equations are conservation laws 

dtp + dif = . (3) 

Constitutive relations express f in terms of p and its derivatives. For instance, if we take 
f = —Dd^p (Fick's law), we get 

dtp - Dd.d'p = . (4) 

Writing p{k, t) = J cf3^g-«A;-a;p(x,t) lioyQ 



p{k,t) = e-^''''p{k,t = 0) . (5) 



This is the characteristic behavior of hydrodynamic modes. It has a hfe T{k) = jyip, which is 
infinite in the long wavelength limit A; ^ 0. Indeed, since p is a conserved quantity it cannot 
disappear locally but can only relax slowly over the entire system. 

The dispersion relation for the hydrodynamic mode is a; = —iDk"^. It shows up as a pole in 
the retarded correlation function 

= (p(f,i),p(0,0))eg , (6) 

where {■■■)eq is thermal equilibrium average, and we assumed {p{x,t))cq = 0. S{x,t) describes 
the fluctuations of the charge density p and can be measured experimentally. 

Relativistic Hydrodynamic: We will work in four flat space-time dimensions with a Lorentzian 
metric rjni, = diag{—l, 1, 1, 1). Define the hydrodynamic expansion parameter (Knudsen num- 
ber) Kn as Kn = Icor/L <^ 1 where Icor is the correlation length of the fluid, and L is the 
characteristic scale of variations of the macroscopic fields. The hydrodynamics equations are 
conservation laws 

a^T^- = 0, d^J^ = , (7) 

where T^*^ is the stress-energy tensor and are symmetry currents. The constitutive relations 
express T'^'^ and in terms of the energy density e{x), the pressure p{x), the charge densities 
Pa{x) and the four-velocity field u'^{x) = (7,7/3*) satisfying u^u^^ = —1. 

Consider for simplicity neutral hydrodynamics. The constitutive relation has the form of a 
series for the stress-energy tensor, 

oo 

T'^'^(x) = ^rj;'(x), (8) 

where T^'^j^ ~ [Knf . The expansion is a derivative expansion, where / counts the number of 
derivatives. Keeping only the first term Z = in the series gives ideal hydrodynamics and the 
stress-energy tensor reads 

^J^ = euV+pP''^ (9) 

where P*^^ = r]^^ + u^vy . Since there are four conservation equations and five fields e,p,u^, 
one more equation is needed. This is called the equation of state e{p). For instance, in CFT 
hydrodynamics: = and the equations of state reads e = 3p. In this case, there is only one 
dimensionfull quantity, the temperature T, and e,p ~ T^. The ideal conformal fluid stress-energy 
tensor can be recast in the form 

T^'^^ = r^[,yA''^ + 4u^nn , (10) 

up to and overall constant coefficient. The entropy current is defined as = su^^ where s is 
the entropy density s ^ T'^, and is conserved in ideal hydrodynamics d^^S^ = 0. 

The first order dissipative hydrodynamics is obtained by keeping also the I = 1 term in the 
series. Going out of equilibrium, there is an ambiguity in the definition of the fields, which 
requires a choice of a frame. In the case of neutral hydrodynamics there is one such frame, the 
Landau frame u^T^'^j = 0. In this frame, the local rest frame of the flow is where the energy 
density is at rest. The stress-energy tensor reads 

^m = -W^'-C(5a^^")^'^"' (11) 



where 



(12) 



is the shear tensor. There are two transport coefficient at this order, the shear viscosity r] and 
the bulk viscosity Note, that in the conformal case C = and the first order dissipative hydro- 
dynamics of a CFT is determined by only one transport coefficient. The transport coefficients 
can be calculated form the retarded Green functions of the microscopic thermal field theory 
using linear response theory and the Kubo formula. 

The entropy current is no longer conserved at the first viscous order and we have 

d,S^ = ^a,,cT^' + ^{d,un' > . (13) 
3 The Fluid/Gravity Correspondence 

Consider the five-dimensional Einstein equations with a negative cosmological constant 

Emn = Rmn + ^9mn = 0, i? = -20 . (14) 

These equations have a thermal equilibrium asymptotically AdS solution, the boosted black 
brane. The metric in the Eddington-Finkelstein coordinates reads 

gl^ldy^'dy'' = -lu^dx^'dr - f[br]Uf,u^dx''dx'' + r'^P^^dx^'dx" , (15) 

where y = {x^,r), is the 4-velocity boost vector, T = l/vri is the temperature and /(r) = 
1 — \. The metric (jlSp is a solution of equations ()14p when b and n'^ are constants. It has an event 
horizon (null hypersurface) with planar topology at r = , whose normal is {n^,rf) = (n^,0). 
The Bekenstein-Hawking entropy density s = ^ and the normal combine to give the entropy 
current = -^u^, where we set Gjv = 1- 

The metric (|15p is no longer a solution when 6(x"), M^(a;"), and needs to be corrected. One 
looks for a solution of the Einstein equations ([H]) by the method of variation of constants 

9mn 9mn ~^ 9mn 

+ ... , (16) 

where g'^n includes first derivatives of 6(x") and n'^(x"). Imposing AdS asymptotics and regu- 
larity at the horizon, one finds that the momentum constraint Einstein equations -EJ^ = give 

the boundary CFT hydrodynamics equations d^T^'^ = as a series expansion in derivatives^. 
The Hamiltonian constraint i?^ = is the equation of state. The ratio of the shear viscosity 
to the entropy density that is calculated from the gravitational background takes the universal 
value 2 = ^. 

The CFT hydrodynamics equations can be viewed also as the Gauss-Codazzi equations 
governing the dynamics of the event horizon^^. In particular, the n ull h orizon focusing equation 
is equivalent to the entropy balance law of the hydrodynamic fiuid^^^. 

One can generalize the fiuid/gravity correspondence to charged and non-conformal hydro- 
dynamics by including bulk gauge A'^ and scalar (pi fields, respectively. Now the bulk viscosity 
is nonzero. Remarkably, one can derive a simple formula for the ratio of the bulk to shear 
viscosities^Sl 

c__ 




(17) 



where s is the entropy density, are the charges associated with the gauge fields j4^, are 
the values of the scalar fields on the horizon, and the derivatives are taken with couplings and 
mass parameters held fixed. 



4 Quantum Anomalies 



The hydrodynamic description exhibits interesting effects when a global symmetry current 
of the microscopic theory is anomalous 

d,Jli = lCabce'-"'''F^,F'^^ . (18) 

Cabc is the coefficient of the triangle anomaly of the currents Jq and J^. The form of an 
anomalous symmetry hydrodynamic current is modified by a term proportional to the vorticity 
of the fluid 

u^" ^ le^'^'u^dxu, . (19) 

This has been first discovered in the context of the fluid/gravity correspondence^^. 
The global anomalous symmetry hydrodynamic current takes the form 

.■M _ + (e^ - TP^-d, + Ca^^' + Cif ^ B'f' . (20) 



= PaU- 



Ej^ = F^^Uu, B"-^^ = ^e^^^^^UyF^p, while pa, T, fia and are the charge densities, temperature, 
chemical potentials and the conductivities of the medium. The anomaly coefficients can be 
calcul ated from the requirement that the entropy current has a positive divergence, d^S'^ > 
They read 



b,.c , na rr2 '^Pa ( ^ ,.b,,c,,d , na ..brr2 



e + p \3 

ti^) _ n ,,c f^ry. ...c.d , «.^2 



Cat' = a6c/^'=-^(^2^6cdMV + /36T^j • (21) 

The coefficients /?„ are related to gravitational anomaly that is to the triangle anomaly 
diagram of the current and two stress-energy tensors. 

Consider next possible experimental signatures of the axial current triangle diagram anomaly 
in a hydrodynamic description of high density QCD. In RHIC and LHC the chemical potentials 

are small compared to the energy density, hence = CabcP fJ-'^^Cab — CabcfJ-'^- The chiral 
magnetic effect and the chiral vortical effect correspond to charge separation and baryon number 
separation, respectively^'. 

The chiral magnetic effect corresponds to the generation of an electric current in the direction 
of the magnetic field j haB, using Cabc of two electric currents and one axial current in 
(|20p. with being the axial chemical potential. The chiral vortical effect corresponds to the 
generation of a baryon number current in the direction of the vorticity vector j ~ haPb'^^ using 
Cabc of two Baryon number currents and one axial current in (|20p . with /i^ being the Baryon 
chemical potential. 

Another possible experimental signature is based on the idea that the axial charge density, 
in a locally uniform flow of massless fermions, is a measure of the alignment between the fermion 
spins. When the QCD fluid freezes out and the quarks bind to form hadrons, aligned spins result 
in spin-excited hadrons. The ratio between spin-excited and low spin hadron production and its 
angular distribution may therefore be used as a measurement of the axial charge distribution. 
One predicts an enhancement of spin-excited hadron production along the rotation axis of the 
collision, the cone compared to the belt (see flgure 2)1^. 

New chiral effects arise in superfluid hydrodynamics ^^^^I^, that is when there are sponta- 
neously broken symmetries. This is a relevant framework for QCD at high densities and low 




Figure 2: An enhancement of spin-excited hadron production along the rotation axis of the colhsion. 

temperature^^. One possible observable effect is the chiral electric effectJ^^, which is a generation 
of an electric current perpendicular to the electric field 

J^'ee = cV'^'u^epEl , (22) 

compared to the standard electric conductivity term Jc^onduct ~ ^""^E^- ^abc are the transport 
coefficients proportional to the triangle anomaly diagram coefficient Cabc , is the phase gra- 
dient of the broken symmetry, which is proportional to the velocity of the superfluid part. If 
the broken charge is gauged as in the superconducting case, the gauge fields will be dynamically 
excluded from the bulk of the system. When we have unbroken gauged charges, then their 
associated gauge fields may enter the superfluid. 



5 Turbulence and Singularities 

The Reynolds number TZe is a dimensionless parameter that measures the relative strength of 
the ideal and the viscous parts of the stress-energy tensor, and is inversely proportional to the 
Knudsen number. For instance, in CFT hydrodynamics we get 

TT 

Tie ~ — , (23) 

rj/s 

where T is the temperature, L is the characteristic scale of variation of the macroscopic fields 
and ij/s is the ratio of the shear viscosity and the entropy density. For a large Reynolds number 
TZe ~ 10^ the fluid exhibits turbulence. 

Non-relativistic turbulence exhibits a universal structure as we will see in the next section. Is 
there a universal structure also in relativistic turbulence? Consider the hydrodynamics equation 
with a random force term 

d'T^u = U ■ (24) 

In the inertial range of scales, scales much smaller than the scale of the force and much bigger 
than the viscous scale, one can derive the exact scaling relation^ISI 

{To,iO,t)T,,{r,t)) = eri , (25) 

where (...) is an average with respect to the random force and e is a nonuniversal constant. 
Similarly, in charged hydrodynamics with a conserved symmetry current one can derive 



{Jo{0,t)Ji{r,t)) = en 



(26) 



It would be interesting to look for experimental setup where these can be verified. For gold 
collisions at RHIC, the characteristic scale L is the radius of a gold nucleus L ~ 6 Fermi, the 
temperature is the QCD scale T ~ 200 MeV, and ^ r-j -1- is a, characteristic value of strongly 
coupled gauge theories. With these Re is too small for an experimental realization of relativistic 
turbulence. An experimental setup, where one may be able to study universal properties of 
relativistic turbulence is condensed matter physics. For instance, it was noticed that there is 
an emergent relativistic symmetry of electrons in graphene near its quantum critical point, for 
which a relativistic nearly ideal fluid description may be appropriate. 

An important issue in the hydrodynamic description is whether starting with appropriate 
initial conditions, where the velocity vector field and its derivatives are bounded, can the system 
evolve such that it will exhibit within a finite time a blowup of the derivatives of the vector field. 
Physically, such singularities if present, indicate a breakdown of the effective hydrodynamic 
description at long distances and imply that some new degrees of freedom are required. The 
issue of hydrodynamic singularities has an analogue in gravity. Given an appropriate Cauchy 
data, will the evolving space-time geometry exhibit a naked singularity, i.e. a blowup of curvature 
invariants and the energy density of matter fields at a point not covered by a horizon. For a 
first step in relating the hydrodynamics and gravity singularities using the Penrose inequality 
I2Q], seem 



6 Nonrelativistic Flows 



Consider the equations of relativistic CFT hydrodynamics 

a^n^' + SDlnT = ^a^^a^"^ , 
a„ + P^^^^lnT = ^pi^{d^a^ - ScT^aJ , (27) 

where D = u^d^, = Du^ and we took ^ = i-. 

The incompressible Navier-Stokes equations can be obtained in the nonrelativistic limit of 
relativistic hydrodynamics^^!!. Expand = (1 + + ...,f*),T = To(l + P + ...) where we 
scale ~ e, 5i ~ e, dt ~ e^, P ~ e^, where e ~ 1/c. 

The first equation in (|27p gives the incompressibility condition Sjf * = 0. The second equation 
gives 

dtv' + v^djv' = -d'P + v^v' , (28) 

where u = is the kinematic viscosity. 

Most nonrelativistic fluid flows in nature are turbulent. The Reynolds number takes the 
form 

LV 

ne = — , (29) 

where L and V are, respectively, a characteristic scale and velocity of the flow. TZe is generically 
large since the kinematic viscosity is generically low. For instance, the kinematic viscosity of 
water at room temperature is ~ 

There is experimental and numerical evidence that in the inertial range of distance scales, 
the flows exhibit a universal behavior 

Sn{r) = (((v(x) - v(y)) • V r^" , (30) 

where r = x — y , and the anomalous exponents are measurable real numbers. The 1941 exact 
scaling result of Kolmogorov .^3 = 1 agrees well with the experimental data. The major open 
problem of turbulence is to calculate the anomalous exponents It would be interesting to 
see if gravity can shed light on this calculation^SS. 



Note, that the relativistic scahng relation (j25p reduces in the nonrelativistic Hmit to Kol- 
mogorov relation 

{vj{0,t)vi{r,t)vj{r,t)) = en . (31) 
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